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Abstract 

An adaptive unstructured remeshlng technique Is 
evaluated for integrated fluid-thermal-structural analysis. 
The technique Is combined with the finite element 
method to solve: (1) the Navler-Stokes equations for 
high-speed compressible flow, (2) the energy equation 
for the structural-thermal response, and (3) the quasi- 
static equilibrium equations for the structural response. 
The remeshlng technique and the analysis solution 
procedure are described. The effectiveness of the 
approach is evaluated with two application studies. The 
flow analysis of Mach 8 shock-shock Interference on a 
three-lnch-dfameter cylinder Is used as the first 
application study to demonstrate the capability of the 
remeshing technique and to examine proper remeshlng 
Indicators for the Invlscld and boundary layer regions. 
The applicability of the approach for the thermal and 
structural analyses of the structure is evaluated In the 
second application study of a 0.25-lnch-dlameter 
convectively cooled leading edge subjected to Intense 
aerodynamic healing. Issues associated with remeshlng 
indicators for thermal stress problems are Identified. 

Nomenclature 


x, y coordinate directions 

X,Y principal coordinate directions 

At time step 

AU U n+1 -U n 

P thermal expansion parameter 

8 boundary layer thickness 

e flow total energy, Eq. (2), or strain 

components, Eq. (18) 

X eigenvalues 

4 key parameter for remeshing, Eq. (20) 

p density 

ox. oy, Txy fluid stress components, Eq. (2), or solid 
stress components, Eq. (7) 

Subscripts 
F fluid 

L left element 

n right element 

T thermal 

s structural 

1 element internal flux 

2 flux across element boundary 


A coolant passage area. Eq. (5), or finite 

element area, Eq. (10) 

[A*] Jacobian matrix, Eq. (8) 

(BJ viscous flux matrix, Eq. (11) 

c specific heat, Eq. (4), or material elastic 

constants. Eq. (18) 

E, F x and y flux components 

H load vector, Eq. (3) 

h convective heat transfer coefficient, Eq. (5), 

or element spacing parameter, Eq. (22) 

[ I ] Identity matrix 

(K) stiffness matrix 

k thermal conductivity 


1, m 

components of unit normal vector 

rh 

coolant mass flow rate 

[M] 

mass matrix 

IN) 

element Interpolation or weighting function 

n 

unit normal vector 

P 

coolant passage perimeter 

P 

pressure 

q 

heat flux 

{R> 

load vector 

s 

edge length of finite element 

T 

temperature 

To 

reference temperature for zero stress 

t 

time 

U 

conservation variable 

u. v 

flow velocity components, Eq. (2), or 
displacement components, tq. (19) 
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Superscripts 

l invlscld 

V viscous 

T transpose 

l,| summation Indices 

n time step index, t » nAt 


Introduction 

Design of lightweight structures and thermal 
protection systems for hypersonic vehicles depends on 
accurate prediction of aerothermal loads, structural 
temperatures and their gradients, and structural 
deformations and stresses. Adaptive unstructured 
remeshlng techniques combined with the finite element 
method have been shown to Improve significantly the 
efficiency and accuracy of the flow analysis. These 
benefits should accrue for both the thermal and structural 
analyses of structures. In addition, a universal 
remeshlng procedure Is needed for the development of 
an integrated fluid-thermal-structural analysis capability. 

Research Is currently underway at NASA Langley 
Research Center to develop an Integrated analysis 
approach for accurately predicting each disciplinary 
behavior and their Interactions. The approach uses the 
finite element method to solve: (1) the Navler-Stokes 
equations for the aerodynamic flow field and the 
aerothermal loads; (2) the structural energy equation for 
heat transfer and temperature distribution; and (3) the 
structural equilibrium equations for deformation and 
stresses. The capability of the approach was 
demonstrated for an undisturbed Mach 6 flow over a 



three-inch-diameter stainless steel cylinder 1 for which 
the experimental dala^ were available for comparison. 
The flow field Is characterized by the bow shock that 
stands-off the cylinder and the thin boundary layer at the 
cylinder surface. Sharp gradients In the flow variables 
occur in these regions and closely spaced elements 
were generated manually to model these flow features. 
The flow solution compares well with the experimental 
data demonstrating the effectiveness of the finite element 
method used. 

For a more complex problem, where a priori 
knowledge of the How physics does not ex!st t an 
appropriate finite element mesh may not be constructed 
easily. An adaptive mesh generation technique is 
required not only to construct and adapt the mesh 
automatically to represent the How field but also to 
reduce the number of unknowns and the analysis 
computational time. Thus the current focus on the 
Improvement of the integrated analysis approach Is to 
develop efficient adaptive unstructured mesh generation 
techniques for each analysis discipline considered 
including their interface requirements. 

Adaptive mesh generation techniques may be 
classified Into two major categories: (1) 

refinement/derefinement, and (2) remeshlng. The first 
category, the adaptive refinement/derefinement 
technique, can be further classified into three 
subcategories: (a) the h method, (b) the p method, and 
(c) the r method. In the h method, the elements in the 
Initial mesh are refined into smaller elements or 
dereffned Into larger elements®, The p method maintains 
the geometry of the elements of the initial mesh but 
increases (or decreases) the order of the polynomial 
used for the element interpolation function*. The r 
method keeps the number of elements and their 
connectivities the same but relocates the nodes®. These 
methods help Increase the analysis solution accuracy but 
each has limitations. Even though nodes could be 
added or removed by the h method, the orientation of the 
original and new elements does not change. Orientation 
of element sides along principal flow gradients Increases 
solution accuracy and reduces the number of elements®. 
Even though the r method can move the nodes for better 
alignment, the method may generate highly distorted 
elements If the number of elements Is fixed. The p 
method may alleviate these restrictions because of the 
flexibility In selecting the order of polynomials on the 
element sides. However, the use of such higher order 
polynomials In combination with the hierarchical concept 
leads to complexities In the formulation of finite element 
equations and Implementation In computer programs. 

The above considerations led to the development 
of the second adaptive mesh category, the remeshing 
technique?. The technique generates an entirely new 
mesh based on the solution obtained from the earlier 
mesh. The technique combined with the finite element 
method has been applied successfully to several 
compressible flow problems with complex flow 
behavior®*®. The purpose of this paper Is to extend the 
adaptive remeshlng technique to both the thermal and 
structural analyses of structures. The adaptive 
remeshlng technique will be evaluated to assess its 
effectiveness for the integrated fluidthermal-structural 
analysis. The governing equations for the aerodynamic 
flow, structural heat transfer, and structural response will 
be Introduced first. The finite element solution algorithms 
for solving these equations will be described. The 


basic concepts of the adaptive remeshlng technique will 
then be explained. Selection of Indicators used for 
construction of new meshes for the three analysis 
disciplines will be discussed. Finally the adaptive 
remeshlng technique will be evaluated for: (1) the fluid 
analysis of Mach 8 shock-shock Interference on a three- 
Inch-dlameter cylinder where the experimental data Is 
available, and (2) the thermal-structural analysis of a 
0.25-lnch-dlameter convectlvely cooled leading edge 
subjected to Intense aerodynamic heating simulating 
Mach 16 flight condition. The Issues of mesh continuity 
along the fluid/structure Interface and between thermal 
and structural analyses for direct Interdisciplinary data 
transfer will also be discussed. 


Integrated Fluld-Therm8l-Structural Procedure 
Governing Equations 

Aerodynamic Flow. The equations for a laminar 
compressible flow are governed by the conservation of 
mass, momentum, and energy. These equations are 
written In conservation form as 

* If - 0 (1 > 


where the subscript F denotes the fluid analysis. The 
conservation variables vector, {Up), and the flux vectors 
In the x and y directions, (E F ) and {Fp}, are given by 
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Superscripts l and v represent the Invlscid and viscous 
flux vector components, respectively. The pressure p Is 
related to the total energy assuming a calorically perfect 
gas (constant ratio of specific heats). The stresses a,, 
o y , and t*y are related to the velocity gradients by 
Stokes' hypothesis. The heat fluxes q x and q y are 
related to the temperature gradients by Fourier's law. 
The temperature-dependent viscosity Is computed from 
Sutherland's law and the thermal conductivity Is 
computed assuming a constant Prandtl number of 0.72. 

Structural Heat Transfer. The thermal response of 
the structure Is governed by the energy equation which 
1$ written In conservation form as 

S< u t>*R< e t>*£< f t> - h t < 3 > 

where the subscript T denotes thermal analysis. For 
transient heat conduction, the conservation variable, Ur, 
and the heat flux components, Et and Ft, are 
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(8) 


dU T =. p s c # dT, 

Et - 0. - ".fir 

f t * <V ■ A 37 « 

and Ht Is the heat source. The structural heat flux 
components q* and qy are related to the temperature 
gradients by Fourier's law. 

For a structure with Internal convection cooling, 
the energy equation for the coolant flow In the local x- 
direction based on assumed uniform cross-sectional or 
bulk temperature 1 can also be written In the 
conservation form of Eq. (3) where 

dU T - p F c F dT F 

E y - mCpTp/Ap - k p 3T/9x 

F y - 0 

H t = hP(T,-T F ) (5) 

The flux Et consists of the energy transport by 
convection (1st term) and conduction (2nd term), and the 
flux Hr represents the heat transfer between the 
structure and Ihe coolant. 


^ 4 E R ♦ I A* I (U L - U R )1 


where |A*| denotes the absolute Jacobian matrix 
associated with the element flow velocity components 
and the local speed of sound. Details of the Jacobian 
matrix A* are given In Ref. (12]. 

In the computation of viscous flux components 
such as the heat flux In the x-dlrectlon (q« = -kF 3T/3x), 
the element temperature gradient is computed from 


^ 7 ^ - C N <x.y) 1 I U > 0) 

where (N(x,y)J Is a matrix of finite element inlerpolatlon 
functions. The vector on the right hand side of the 
equation contains the nodal temperature gradients, 
determined at node K, by 

fr'K - 5T-( J BT N K ' |S - «»> 

K * A 


where n is the unit normal vector along the element 
Interface s and M« Is the lumped mass at node x. 

Increments In time of the element conservation 
variables, All - U n+1 - U n , are determined from 


Structural Response. The structural response Is 
governed by the quasi-static equilibrium equations 
assuming that the inertia effect is negligible. These 
equations are written in conservation form as 
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The flux vector components, (E s ) and (F S J, are 

{E.) T - (o, t^] 

(f/ - I V 0 y ] (7) 

where the stress components, o*, cr y , and x« y are related 
to the displacement gradients and the temperature by 
generalized Hooke's law. 

Finite Element Solution Algorithms 


The matrix [B], which Is associated with the viscous 
fluxes Is given In Ref. (10], accelerates the solution 
convergence In the viscous dominated regions. This 
fluid algorithm provides first order of accuracy (0(h)) In 
the Inviscid flow dominated region but the order of 
accuracy Is close to 1.5 (0(hi *)) In the viscous flow 

dominated region 1 9. 

Flux-Based Taylor-Galerkln Algorithm. The basic 
concept of the flux-based Taylor-Galerkln algorithm Is to 
express the element flux components E a nd F In the 
form 


E(x.y.t) 

F(x.y.t) 


[N(x,y)l 

(N(x.y)( 


( E ( t )} 
<F(t)) 


( 12 ) 


An Impllclt/expllclt upwind cell-centered algo- 
rithm®, 10 is used to solve the Navier-Stokes equations, 
Eqs. (1)-(2). A flux-based Taylor-Galerkln finite element 
algorithm®. 1 1 Is used to solve the thermal and structural 
equations, Eqs. (3)-(7). For brevity, only essential 
features of the algorithms are highlighted herein. 

Upwind Fluid Algorithm. The basic concept behind 
the upwind cell-centered algorithm Is to determine the 
flux across element interfaces by Roe's averaging 

_t 

procedure 1 2. The average Inviscid flux E (the 
"overbar* denotes an average quantity normal to the 
element Interface between the left element, L, and the 
right element, R) Is given by 


where (E), and (F) are vectors of the element nodal flux 
quantities. 

For the thermal analysis, the final finite element 
equations are In the form 

(M](AU T r' - (R t ]" + (R t )J (13) 

where at time n+1, AU n+1 = U n+ i - U n . The mass 
matrix, (M], Is determined from 

(M) « | (N) (N) dA (14) 

A 

The two vectors on the right hand side of Eq. (13) are 
associated with the fluxes within the element (subscript 
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1) and across the element boundary (subscript 2). 
These vectors are defined by 

{R t )J = At (J{|l}{N]dA {E t }" 

A 

+ J(M,( N ] d A (F T } n ) (15) 

A 

{ R t >2 = -At J(N)(N)ds (l|E T )% m|F T } n ) 

8 

= -At J(N)[N)ds (q)" (16) 

s 

The vectors (Et) and {F T } contain nodal heat fluxes 
which depend on the nodal temperature gradients. The 
element nodal temperature gradients 3T/3x and 3T/dy t 
and all the finite element matrices (Eqs. (13) (16)) can 
be evaluated In closed form (Le., numerical Integration is 
not required) for any type of finite elements Including the 
4-node quadrilateral and the 8 node hexahedral 
elements. 

For the structural analysis, the final finite element 
equations are In the same form as the thermal equation, 
Eq. (13), but without the transient term, 

* ( R, ) 2 - 0 (17) 


model large solution gradients. Based on these ideas, 
the adaptive remeshing technique consists of two steps: 
(1) the determination of the new element sizes and their 
orientations, and (2) the construction of the new mesh. 

In the determination of the element size and 
orientation, the solid mechanics concept of determining 
the principal stresses and their directions from a given 
state of stress at a point is employed. For example In the 
thermal analysis, small and clustered elements are 
needed in regions of large temperature gradients. The 
temperature is thus considered as a key parameter to 
indicate where clustered elements are needed. At a 
typical node In the old mesh, the second derivative 10 of 
the key parameter with respect to the global coordinates 
x and y are first computed, 



where $ denotes the key parameter used for remeshlng. 
Then the eigenvalues which represent the principal 
quantities In the principal directions X and V where the 
cross derivatives vanish are determined, 


where these vectors are Identical to Eqs. ( 1 5)-( 1 6) 
except that the subscript T is replaced by the subscript s 
everywhere. The element nodai flux vectors {E$} and 
(F s ) (similar to (Et) and (F T J in Eq. (15)) now consist of 
the stress components a x , a Y and t xy . By using the 
thermal stress-strain constitutive relations, 

o, * + P, (T-T 0 ) IJ- 1.2,3 (18) 

and the strain-displacement relations, these nodal stress 
components associated with the (R # }i vector can be 
written in terms of the nodal displacement components. 
Further algebraic manipulation of the finite element 
equation, Eq. (17), results In the standard form of 

IK)(U} = {Rj j 2 + { r t , (19) 

where (U) consists of the unknowns of the nodal 
displacement components. [KJ Is the element stiffness 
matrix, and (R-r) Is the thermal load vector. These 
matrices are evaluated in closed form and the details 
are given In Ref. [11]. 


Adaptive Remeshlng Procedure 
Remeshlng Concept 

The basic Idea of adaptive remeshlng Is to 
construct a completely new mesh based on the solution 
obtained from the previous mesh. The new mesh will 
consist of clustered elements In regions with large 
solution gradients and few elements In the regions where 
the gradients are small. Element orientations are In 
principal directions to provide the most accurate solution 
with a minimum number of elements. As an example, the 
shorter element side Is In the direction normal to the 
shock line or through the boundary layer thickness to 



These eigenvalues are the remeshlng indicators and are 
used to determine the element spaclngs h, and h 2 In the 
two principal directions using the condition of 


hjX, = h^ - constant (22) 


which has been shown to distribute the solution error 
equally, a condition tor an optimal mesh 13. This process 
Is performed for all the nodes In the old mesh leading to 
the final condition 7 of 


h*X 1 m h*X 2 * constant *■ h mh Km < 23 > 


Using this condition and specifying the required 
minimum element spacing hmio, the new element 
spaclngs based on the solution of the old mesh are 
obtained and the new mesh is constructed. 

Mesh construction is based on an advancing front 
technique 7 . Nodes are first placed along the domain 
boundary. Spaclngs between these nodes are dictated 
by the element spacing parameter h obtained earlier. At 
this Instance, the front consists of a sequence of straight 
line segments which connect consecutive boundary 
nodes, l.e., the domain boundary. Elements adjacent to 
the domain boundary are then constructed. The sizes 
and orientations of these elements are guided by the 
computed element spacing parameter h and the principal 
directions. As new elements are constructed from the 
domain boundary, the front Is updated and represented 
by the new element sides. As the mesh construction 
goes on, the front changes Its shape continuously and 
propagates toward the domain Interior, The generation 
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process Is terminated when the domain Is tilled 
completely with elements and the front vanishes. 

Based on the condition shown_in_Eq. (23), the 
element size Is generated according to the given 
minimum element spacing h m i„. Specifying too small a 
spacing h m i n could result In a model with a large number 
of elements which may be impractical to perform the 
analysis. On the other hand, specifying too large a 
spacing h m i n may result In an Inadequate analysis 
solution accuracy, or excessive analysis and remeshlng 
cycles to achieve the required solution accuracy. These 
factors must be considered prior to generating a new 
mesh. Note that, because the technique generates an 
entirely new mesh with different nodal locations from the 
old mesh. Interpolation of the solution from the old to the 
new mesh is necessary. The Interpolated nodal 
quantities are used as the initial conditions for the new 
mesh to increase the analysis solution convergence. 

Flufd-Thermal-Structural Remeshlng Parameters 

The adaptive remeshlng technique described 
requires a selection of proper key parameters (<|> In Eq. 
(20)) for remeshlng so that generaled elements are 
clustered where needed. Selection of the key 
parameters depends on the analysis discipline and its 
applications. For thermal problems, It Is obvious that the 
temperature should be used as the key parameter so that 
the mesh generated can represent steep temperature 
gradients. For structural problems, stress Is an 
appropriate choice for Ihe key parameter so that regions 
with high stress concentration will be represented. 
However, the key parameter representing the stress 
should be a scalar quantity such as the Von Mises stress 
defined by, 

CT VonM.*» “ j= Y K-V 2+ + + 6t2 *Y < 24 > 

Note that the key parameter selected for remeshlng may 
vary with applications. For problems that require 
accurate deformations rather than the stresses, an 
absolute displacement quantity may be preferred as the 
key parameter. 

For high-speed compressible flow problems, 
selection of Ihe proper key parameters for remeshlng 
may be more difficult. In the Inviscld flow field, the fluid 
density Is normally selected as the key parameter^ 
because density exhibits high gradients across shock 
and flow expansion waves. In the viscous dominated 
flow regions, such as the boundary layer, other key 
parameters may be more appropriate. Key parameters of 
heating rate and skin friction may be used for more 
accurate analysis prediction of the aerolhermal loads. 
The selection of these key parameters for remeshlng In 
the three different analysis disciplines will be discussed 
in detail in the next section. 


Evaluation of Integrated Analysis and 
Adaptive Remeshlng Procedures 

Two applications are presented to assess the 
adaptive unstructured remeshlng procedures for 
integrated fluid-thermal-structural finite element analysis. 
The fluid analysis of Mach 8 shock-shock Interference on 
a three-inch-diameter cylinder Is used as the first 
application to Illustrate the successful Implementation of 


the adaptive remeshlng technique for complex flow and 
to Investigate the meshing issues for the Integrated fluid- 
thermal-structural model. The applicability of the 
remeshlng technique for thermal and structural problems 
Is demonstrated In the second application for a 0.25- 
Inch-diameter convectively cooled leading edge 
subjected to Intense aerodynamic heating. Remeshlng 
parameters for thermal and mechanical stress 
predictions are Identified. 

Mach 8 Shock-Shock Interference On a Cylinder 

Leading edges of hypersonic vehicles that 
experience Intense stagnation point pressures and 
heating rates are a significant challenge to the designer. 
For engine leading edges, such as the cowl shown In 
Fig. 1 , intense aerothermal loads occur when the cowl 
bow shock Is Intersected by an oblique shock resulting in 
a supersonic jet that Impinges on the leading edge 
surface. The experimental configuration (lower left of Fig. 
1), which simulates the vehicle forebody and cowl 
leading edge, was used to define the aerothermal 
loads* 14. The experimental configuration Is rotated 
180° relative to the vehicle. The schlleren photograph 
shows the supersonic jet interference pattern Impinging 
on the surface of the cylinder. The Interference pattern 
produces Intense local amplification of the pressure and 
heat transfer rate In the vicinity of the jet Impingement. 
The undisturbed (absence of incident oblique shock and 
Interference pattern) stagnation pressure and heating 
rate can be amplified by factors from 6 to 30 depending 
on the shock strength and the free stream Mach 

number2,14. 

The flow Interaction pattern superimposed on the 
flow computational domain and the three-lnch-dlameter 
cylinder is shown schematically In Fig. 2. The Inflow 
conditions above and below the oblique shock are: (1) 
Mach 8.03 flow at an angle of attack of zero degrees 
(a=0°) and static temperature of 200 "R, and (2) Mach 
5.25 flow at an angle of attack of 12.5 degrees (a=12.5°) 
and a static temperature of 430 °R. The supersonic jet 
Impinges on the cylinder surface approximately 22 
degrees below the cylinder horizontal center line. 

A first mesh shown In Fig. 3 Is generated using a 
background mesh concept described In Ref. (7). A total of 
4192 triangles are generated In the Inviscld region and 
1568 quadrilateral elements In the boundary layer. The 
reasons for using quadrilateral elements in the boundary 
layer will be described In detail below. Using this mesh, 
the fluid analysis Is performed to obtain a flow solution as 
illustrated by the Mach number contours shown In Fig. 4. 
Based on this flow solution and the use of the fluid 
density as the key parameter for remeshlng, the second 
mesh Is created as shown In Fig. 3. The flow solution 
obtained from the first mesh Is Ihen Interpolated and 
used as the Initial condition for the second mesh. The 
same procedure Is repeated on the subsequent meshes 
until the converged flow solution Is achieved (a total of 
three meshes In this case). 

The Mach number contours shown In Fig. 4 
demonstrate the Improvement of the solution quality with 
the adapted meshes. The Mach number distribution 
obtained from the third mesh shows Improved sharpness 
of the shock Interference pattern. The contours show a 
supersonic jet (Mach number - 2) submerged within the 
subsonic regions between the bow shock and thp 
cylinder. The supersonic flow In the jet terminates 
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through a nearly normal shock prior to Impinging on the 
cylinder surface. 

The analytically predicted surface pressure 
distribution from the third mesh Is compared with the 
experimentally measured pressures 2 In Fig. 5. The 
predicted and experimental pressures are normalized by 
the undisturbed stagnation pressure (po = 10.61 psla). 
Good agreement between the analytical and 
experimental results are obtained for the pressure 
distribution, the peak pressure and its locations. Details 
of the finite element mesh and the flow temperature In the 
interaction region are shown in Fig. 6. On both sides of 
the supersonic jet, the fluid temperature increases 
abruptly across the bow shocks from a relatively low 
temperature (200 °R and 430 °R) to approximately 2,700 
°R. The temperature gradients in the shock layer (region 
between the bow shock and the cylinder) are relatively 
small except In the thin boundary layer where the 
temperature drops sharply to the cylinder surface 
temperature of 530 °R. Inside the supersonic Jet, the fluid 
temperature Increases slightly from the free stream 
temperature to approximately 1,000 °R. As the Jet stream 
approaches the cylinder surface, the fluid temperature 
Increases abruptly across the jet normal shock to 
approximately 3,000 fl R in a small stagnation region next 
to the cylinder surface and drops sharply to the cylinder 
temperature of 530 6 R resulting In high aerodynamic 
heating rates. 

The major difficulty In the fluid analysis Is the 
prediction of the aerodynamic heating rates because a 
very accurate resolution of the flow temperature gradient 
normal to the cylinder surface Is required. For the 
undisturbed Mach 8 flow, the predicted stagnation 
heating rate of approximately 42 Btu/ft 2 -sec, which was 
obtained from a finite element solution 1 and a viscous 
shock layer solutionis, is lower than the experimental 
value of 61 .7 Btu/ft 2 -sec (Ref. [2]). The predicted and the 
experimental interference heating rates are, therefore, 
normalized by their respective undisturbed stagnation 
olnt heating rates and are compared In Fig. 7. The 
eating rate distributions are similar and the peak 
heating rate locations agree. However the predicted 
peak heating rate amplification Is about 60% of that 
measured from the experiment. This difference may be 
attributed to several sources including the finite element 
mesh, the analysis algorithm, and the mathematical 
model. The finite element mesh near the boundary layer 
may be Inadequate to represent the steep temperature 
gradients accurately. The flow analysis algorithm, Eqs. 
(8)-(11), provides an accuracy on the order of the 
element characteristic length h raised to the 1.5 power 
(0(h 15 )). In addition, the fluid flow mathematical model 
does not account for any turbulence. As discussed In 
Ref. 2, the higher experimental heating rate may be 
attributed to free stream turbulence In the test stream 
which emanates from (1) the turbulent boundary layer on 
the nozzle and/or (2) shear layers that separate the 
supersonic Jet from trie subsonic regions. 

As mentioned earlier, layers of quadrilateral 
elements are used In the boundary layer. Detail of these 
quadrilateral elements are shown in Fig. 8. The use of a 
structured mesh of quadrilateral elements in the 
boundary layer provides several advantages In the flow 
analysis. One advantage Is that slender (high aspect 
ratio) elements can be used with thefr longer sides in the 
flow direction where the solution gradients are small and 
their shorter sides normal to the surface to model steep 
solution gradients through the boundary layer thickness. 


A second advantage Is that the artificial dissipation 
required for the numerical stability, which could smear 
the solution gradients, can be controlled easily or turned 
off completely 1 ®. Although the boundary layer thickness 
varies along the cylinder, the thickness (5 = 0.06 Inches 
In Fig. 8) is assumed constant for the purpose of 
constructing the boundary layer mesh. Across the 
boundary layer thickness, 32 layers of graduated 
quadrilateral elements are used with a geometric 
stretching factor of 1.25. The thickness of the first 
element at the surface Is l.xlO 5 Inches and was chosen 
such that It should model the temperature gradient 
expected from the experimental heating rate 2 . 

For general problems, the boundary layer 
thickness Is not known and varies In the flow direction. 
An approach to estimate the boundary layer thickness 
using the fluid heat transfer as the key parameter for 
remeshlng has been proposed recently 10 . The 
approach Is Illustrated In Fig. 9. The boundary layer 
thickness, 5, Is estimated from the absolute quantity of the 
heat flux vector, | q |, which Is a maximum at the surface 
and Is Insignificant above the edge of the boundary layer. 
The size and distribution of the quadrilateral element 
layers are guided by the magnitude of the flow shearing 
stress (t,(y is proportional to 9u/c)y) which Is a maximum at 
the cylinder surface and decreases quadratically through 
the boundary layer thickness. Based on this Idea, an 
adaptive remeshlng technique using quadrilateral 
elements In the boundary layer can be constructed and 
combined with the current capability of the adaptive 
remeshlng technique for triangular elements in the 
invlscid flow region. 

To generate an integrated fiuid-thermal-structural 
finite element model, the use of common nodes along the 
fluid/structure Interface Is preferred 1 ® as highlighted In 
Fig. 10. This approach not only permits direct data 
transfer between the different analysis disciplines but 
also provides consistency for the analysis formulation. 
As an example of a coupled fluid/thermal problem, the 
interface nodal temperatures can be obtained by solving 
the coupled energy equation of the flow and the 
structure 1 with the requirements that at the Interface: (1) 
the temperatures of the fluid and the structure are 
Identical, and (2) the heat flux Is continuous. 

For some problems, the requirement of common 
fluid/structure Interface nodes to preserve the mesh 
continuity may result in an excessive number of nodes in 
one of the disciplines. For example, a high stress 
concentration may occur near the fluid/structure Interface 
which may require a large number of the structural 
nodes even though the flow near that region Is simple. 
Several adaptive meshing techniques are feasible for the 
thermal-structural analysis of the structure as Indicated in 
Fig. 10 which include the mesh refinement/derennement 
technique and the remeshlng technique developed for 
the flow problems. Before selecting one of these 
techniques, their effectiveness for thermal and structural 
problems must be evaluated. An Initial evaluation of the 
adaptive remeshlng technique for both the thermal and 
structural analyses Is performed In the next section. 

Convectlvely Cooled Leading Edge 

A 0.25-lnch-diameter convectlvely cooled leading 
edge subjected to intense aerodynamic heating Is used 
to evaluate the adaptive remeshlng technique for both 
the thermal and structural analyses, and to Identify 
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requirements for the remeshlng parameters for thermal 
stress problems. The example problem represents a 
-hypersonic vehicle operating at Mach 16 which causes 
the vehicle nose bow shock to Intersect with the cowl 
leading edge bow shock. The Intersection produces the 
type IV supersonic jet interference pattern similar to that 
described In the first application. Because of the smaller 
leading edge diameter and the higher Mach number 
conditions, the heating rate Is very Intense with a peak of 
nearly 30,000 Btu/ft 2 -sec (see Ref. (1)). 

The leading edge geometry and boundary 
conditions are shown schematically In Fig. 11. The outer 
surface Is subjected to the aerodynamic heating and 
emits radiant energy to space. The Inner surface Is 
convectlvely cooled by the direct Impingement of a sonic 
hydrogen jet stream with an Inlet temperature of 50 ®R. 
The 1350 quadrilateral element mesh shown In Fig. 
12(a) is typical of a structured mesh. The mesh has been 
used previously 1 for predicting transient thermal- 
structural response of the leading edge as the vehicle 
nose bow shock sweeps across during the vehicle 
acceleration. The mesh is graded In the radial direction 
but Is uniform In the circumferential direction. The 
leading edge material Is assumed to be copper because 
of Its high thermal conductivity. The copper material 
properties used In the analysis are highly temperature 
dependent and can be found in Refs. (17-18). The 
predicted steady-state leading edge temperature 
contours are shown In Fig. 12(b). The peak temperature 
of 766 °R is at the jet impingement location where the 
heating rate Is a maximum. In spite of the high thermal 
conductivity of the material, the temperature gradients 
are also high in this region. The problem is, therefore, 
suitable to be used for evaluating the adaptive 
remeshing technique for the thermal analysis of the 
structure. 

The same adaptive remeshing procedure 
described for the flow problem Is applied for the thermal 
analysis of the leading edge. A fairly uniform mesh 
consisting of 885 triangles, shown In Fig. 13(a), is first 
constructed and the thermal analysis is performed to 
obtain the temperature distribution. Then the adaptive 
mesh shown In Fig. 13(b) Is generated using the 
temperature as the key parameter for remeshlng to 
provide clustered elements In the region with high 
temperature gradients. The temperature distribution 
obtained from the second mesh Is almost Identical (a 
difference of 0.2%) to that obtained from the nonadaptive 
mesh shown In Fig. 12(b). However, the number of 
unknowns required for the adaptive unstructured mesh is 
reduce to about 25% of those used In the structured 
mesh. Note also that several eleme nts are required 
through the thickness where the peak heating occurs and 
only one or two elements are required near the end of 
the leading edge. Hence, only an unstructured mesh 
technique would optimize the mesh. 

To evaluate the effectiveness of the adaptive 
remeshlng technique to represent both the thermal and 
mechanical stress fields, structural boundary constraints 
at the upper and lower leading edge sections are 
Introduced. The inner surface nodes at Jhe constraint are 
fixed (u - v « 0, see Fig. 14(a)) to simulate the constraints 
caused by the Internal fins that support the leading edge 
from the inner primary structure. To highlight the 
remeshlng capability for generating clustered elements 
In the regions needed, the rest of the nodes on both 
sections are constrained differently. These nodes are 
free to move on the upper section, but are constrained In 


the horizontal direction on the lower section. In addition 
to the thermal loads, the leading edge is also subjected 
to mechanical loads from: (1) the aerodynamic pressure 
which has a distribution similar to that shown In Fig. 5 but 
with a peak pressure of 1 ,000 psla, and (2) the uniform 
Internal coolant pressure of 1 ,000 psla. 

The adaptive thermal mesh shown in Fig. 13(b) is 
used as the initial mesh for the structural analysis. The 
same analysis procedure is applied by first obtaining the 
structural analysis solution and then remeshlng In the 
remeshlng process, nodal temperatures of the new mesh 
are Interpolated from the nodal temperatures of the 
thermal mesh shown In Fig. 13(b). For thermal stress 
problems, high stresses normally occur In regions of high 
temperature gradients. High stress concentrations may 
also occur at corners or supports even though the 
temperature Is uniform. Thus, to represent both the 
thermal stress and the stress concentration, key 
parameters for remeshlng based on the temperature and 
the Von Mises stress are used. A new adaptive mesh 
generated by using these two remeshlng parameters Is 
shown In Fig. 14(a). Elements are clustered near the jet 
Impingement location for modeling high thermal stresses 
from the high temperature gradients, and at the region 
near the upper support where the stress concentration 
occurs. The corresponding structural response, based 
on a linear elastic behavior assumption, Is shown by the 
circumferential stress contours superimposed on the 
deformed leading edge In Fig. 14(b). The figure shows 
that a peak compressive stress of 25 ksl occurs at the jet 
Impingement location and a high stress concentration 
occurs near the singular point at the upper support. 

To support the Idea that the two remeshing 
parameters should be used simultaneously for 
generating an adaptive mesh for thermal stress 
problems, the same analysis approach Is repeated but 
without using the temperature as the remeshlng 
parameter. The adaptive mesh generated using the Von 
Mises stress alone as the remeshlng parameter is shown 
In Fig. 15(a). Clustered elements are generated In the 
region near the upper support similar to the previous 
case. The mesh near the jet Impingement location Is 
much coarser than the mesh for the previous case and 
does not resemble the mesh for the temperature 
distribution. The corresponding predicted structural 
response of the circumferential stress contours Is shown 
In Fig. 15(b). The contours are similar to those obtained 
In the previous case but are degraded due to the coarser 
mesh In the region of the jet Impingement location. The 
difference between the peak compressive stresses of 
these two cases at the jet Impingement location is 
approximately 15%. Note that for thermal stress 
problems in which the structure Is free to expand, the 
thermal stress distribution normally resembles the 
distribution of the temperature. For such cases, the use 
of the Von Mises stress alone as the remeshlng 
parameter may be adequate. 

The adaptive remeshlng technique for the 
structural analysis is further evaluated for the case in 
which the thermal load is not present. The leading edge 
Is subjected only to the mechanical load which consists 
of the external aerodynamic pressure and the internal 
coolant pressure. The Von Mises stress is used as the 
remeshlng parameter and the adaptive mesh generated 
Is shown In Fig. 16(a). Similar to the previous case, 
clustered elements are In the region near the upper 
support where the stress concentration occurs. The 
corresponding predicted structural response of the 
circumferential stress contours Is shown In Fig. 16(b). 
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The predicted compressive stress of 5 ksl at the Jet 
Impingement location Is much lower than the stress that 
occurred for the case with the thermal load. Obviously, 
superposition of the results for separate mechanical and 
thermal stress analyses would require Interpolation of the 
results. Hence, the adaptive remeshing and analysis 
should be used for the combined load case. 

For coupled fluid-thermal-structural analysis, both 
the thermal and structural response of the leading edge 
affect the flow field. The flow field has to be updated to 
Include the effects of: (1) the change of the leading edge 
surface temperature, and (2) the leading edge 
deformation. The leading edge may deform Into or away 
from the Initial flow field and updating the flow 
computational domain Is necessary. The new flow field 
may govern regions previously occupied by the leading 
edge where the flow Information does not exist. These 
are some of the future Issues which will be encountered 
and have to be clarified prior to applying the adaptive 
remeshing technique to coupled Interdisciplinary 
problems 


Concluding Remarks 

An adaptive unstructured remeshlng technique 
was evaluated for Integrated fluid-thermal-structural 
analysis. The technique generates an entirely new 
mesh based on the solution obtained from a previous 
mesh. The new mesh consists of clustered elements In 
the regions with high solution gradients and few 
elements In the regions where the gradients are small. 
The capability of the remeshlng technique was 
demonstrated for the fluid, thermal, and structural 
analyses. The finite element formulations used for the 
three disciplines were presented. The basic Idea 
behind the remeshlng technique was described and the 
remeshlng Indicators and requirements were Identified. 

Two applications were presented to assess the 
effectiveness of the remeshing technique. The fluid 
analysis of Mach 8 shock-shock Interference on a three- 
inch-diameter cylinder was used as the first application 
to Illustrate the remeshlng technique for complex flow 
and to Investigate the meshing Issues for the IntegraSed 
fluid-thermal-structural analysis. Fluid density was used 
as the remeshlng parameter for generating adaptive 
triangular elements In the Inviscid flow region. In the 
current procedure, layers of quadrilateral elements are 
constructed In the boundary layer to represent accurate 
aerothermal loads. The predicted aerodynamic 
pressure and heating rate are In good to fair agreement 
with the experimental data. An approach for generating 
an adaptive structured mesh of quadrilateral elements In 
the boundary layer was described. The approach can 
be combined with the current adaptive unstructured 
remeshlng technique used in the inviscid region to 
achieve an automated remeshlng capability for the 
entire flow domain. 

The remeshlng technique was extended to the 
thermal and structural analyses of the structure. The 
applicability of the technique to thermal stress problems 
was demonstrated with a 0.25-inch-dlameter 
convectively cooled leading edge subjected to Intense 
aerodynamic heating. Temperature was selected as the 
remeshlng parameter In the thermal analysis so that 
clustered elements are generated In regions of high 
temperature gradients. The technique provides the 
same temperature solution accuracy compared to a 
refined structured mesh but with a significant reduction 


In the number of unknowns. For the structural analysis 
of the leading edge subjected to both thermal and 
pressure loads, remeshlng parameters of the 
temperature and the Von Mises stress are used to 
provide clustered elements In regions of high thermal 
and mechanical stresses as well as regions with a stress 
concentration. 

The results from the applications have 
demonstrated the viability of the adaptive unstructured 
remeshlng technique combined with the finite element 
methods to provide efficient accurate solutions to 
complex flow-thermal-structural behavior. 
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Fla. 2 Oblique and bow shock interaction pattern for 
Mach 8 shock-shock interference on a three- 
lnch-diameter cylinder. 



Fig. 3 Adaptive finite element mesh evolution ol the 
flow domain for Mach 8 shock-shock 
interference on a three-lnch-diameter cylinder. 




Second m«*h 



TWrdnwsh 


Fig. 1 Vehicle schematic, experimental configuration, 
and schlieren of type IV shock-shock 
interference. 


Fig. 4 Flow Mach number contours for the three 

adaptive unstructured finite element meshes for 
Mach 8 shock-shock interference on a three- 
lnch-diameter cylinder. 
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Fia 5 Comparison of surface pressure distributions on Fig. 8 Details of the finite element mesh in the 

a three-inch-diameter cylinder subjected to interaction region and layers of quadrilateral 

Mach 8 shock-shock Interference. elements In the boundary layer. 
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Fig. 6 Details of the finite element mesh and fluid 

temperature contours In the interaction region 
for Mach 8 shock-shock Interference on a three- 
inch-diameter cylinder 



Fig. 9 Typical temperature and absolute heat flux 

vector distributions In the boundary layer and 
adaptive boundary layer mesh with 
quadrilateral elements. 


o E^erlment 



8, degrees 


Fig. 7 Comparison of surface heating rate distributions 
on a three-inch-diameter cylinder subjected to 
Mach 8 shock-shock Interference. 



Fig, 10 Integrated fluid-thermal-structural mesh concop; 
and adaptive meshing options for thermal- 
structural analysis of the structure. 
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(a) Adaptive mesh (b) Circumferential stress 


Fig. 1 1 A schematic thermal-structural analysis model 
of 0.25-lnch-dlameter leading edge with 
boundary conditions. 


Fig. 14 Final adaptive mesh and circumferential stress 
for convectively cooled leading edge subjected 
to thermal and pressure loads. 




(a) Nonadaptive (b) Temperature contours 

structured mesh 


(a) Adaptive mesh (b) Circumferential stress 


Fig. 12 Nonadaptive finite element mesh and 

temperature contours for convectively cooled 
leading edge. 




(a) First mesh (b) Second mesh 

Fig. 13 Adaptive finite element mesh evolution for 

thermal analysis of convectively cooled leading 
edge. 


Fig. 15 Adaptive mesh using Indicator of Von Mises 

stress and circumferential stress for convectively 
cooled leading edge subjected to thermal and 
pressure loads. 



(a) Adaptive mesh (b) Circumferential stress 

Fig. 16 Final adaptive mesh and circumferential stress 
for convectively cooled leading edge subjected 
to pressure load. 
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